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ABSTRACT 

We have performed extensive simulations to explore the possibility of detecting eclipses 
and transits of close, sub-stellar and planetary companions to white dwarfs in WASP 
light curves. Our simulations cover companions ~ 0.3 < R p i < 12 R e and orbital 
periods 2h < P < 15d, equivalent to orbital radii 0.003AU < a < 0.1AU. For Gaus- 
sian random noise WASP is sensitive to transits by companions as small as the Moon 
orbiting a V ~12 white dwarf. For fainter white dwarfs WASP is sensitive to increas- 
ingly larger radius bodies. However, in the presence of correlated noise structure in the 
light curves the sensitivity drops, although Earth-sized companions remain detectable 
in principle even in low S/N data. Mars-sized, and even Mercury-sized bodies yield 
reasonable detection rates in high-quality light curves with little residual noise. We 
searched for eclipses and transit signals in long-term light curves of a sample of 194 
white dwarfs resulting from a cross-correlation of the McCook & Sion catalogue and 
the WASP archive. No evidence for eclipsing or transiting sub-stellar and planetary 
companions was found. We used this non-detection and results from our simulations 
to place tentative upper limits to the frequency of such objects in close orbits at white 
dwarfs. While only weak limits can be placed on the likely frequency of Earth-sized 
or smaller companions, brown dwarfs and gas giants (radius w R jup ) with periods 
< 0.1 — 0.2 days must certainly be rare (< 10%). More stringent constraints likely 
requires significantly larger white dwarf samples, higher observing cadence and con- 
tinuous coverage. The short duration of eclipses and transits of white dwarfs compared 
to the cadence of WASP observations appears to be one of the main factors limiting 
the detection rate in a survey optimised for planetary transits of main sequence stars. 

Key words: method: data analysis - occultations - stars: planetary systems - stars: 
white dwarfs 



1 INTRODUCTION 

In recent years we have witnessed considerable progress in 
the search for extra-solar planets. Since the first detection 
of a 'Hot Jupiter ' arou nd the main-sequence star 51 Peg 
|Mavor fc Quelo3 Il995l ) . the number of extra-solar plan- 
ets has rapidly risen, and currently approaching 500. Most 
of these discoveries are the result of radial velocity (RV) 
searches. More recently, an increasing number of extra-solar 
planets (> 80) have been detecte d by dedicated pla netary 
transit surveys including HATnet (|Bakos et al.ll2004l '), TrES 
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(e.g. iBrown fc Charbonneaul |2000|: iDunham et all 12004 



lAlonso et alj|2004l) OGL E (|Udalski et al.ll2002l , 2003). XO 
jMcCullough et al.ll2005lb and WASP, t he UK Wide- Angle 
Search for Planets ( Pollacco et al.ll2006t ). 



Planet detection via the transit technique involves 
searching for periodic dips in stellar light curves as a planet 
occludes a small fraction of the visible disc of the host star 
once per orbit. Only planets with their orbital planes aligned 
within a few degrees to the line of sight will exhibit a transit, 
the probability of such an alignment being around 10% for 
typical 'hot Jupiter' systems. This introduces a constraint 
on the number of observable systems and explains the rel- 
atively low number of transiting planets when compared to 
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radial velocity studies. Importantly, when combined with RV 
measurements, planetary transits offer the unique possibility 
of deriving both the planet mass and r adius, since fo r these 
systems the inclination i is well-known jSacketdll999h . For a 
given planetary radius, the transit depth is directly propor- 
tional to (Rp/R*) 2 , where R p , and R* are the planetary and 
stellar radii respectively. Therefore, planets orbiting solar- 
type stars have extremely shallow eclipses, blocking ~1% of 
the light for a giant planet and ~0.01% of the light for an 
Earth-sized planet. Current ground-based wide-field surveys 
can achieve the necessary photometric accuracy of better 
than 1%, only for the brightest stars (V ~ 9-12 in the case 
of WASP), so the bulk of the planets discovered by transit 
surveys around main-sequence stars have radii in the range 
7?p ~ 0.9-1.8 Rjup. To date the smallest extra-solar planet 
detected in a ground-based transit survey is HAT-P-llb, a 
Nept une-size planet (R v = 0.452 R jup ) transiting a K dwarf 
star (|Bakos et alj|20ld ). 

A major advantage over main sequence primaries is of- 
fered by white dwarf stars. White dwarfs (WDs) are compact 
degenerate objects with R W d~ lR®(Earth radius), and 
represent the final stage of evolution of main— sequence stars 
with masses ^ 8M@(i.e. ~ 97% of all stars in our galaxy). 
Any sub-stellar or gas giant companion orbiting the star, will 
completely eclipse it, while bodies as small as the Moon will 
display relatively large transit depths (~ 3%), with the only 
caveat being that it remains unclear as to whether any such 
systems survive beyond the latter stages of stellar evolution. 
The strong gain in the planet— to— star relative dimensions 
opens up the possibility of detecting low— mass sub— stellar 
and in particular terrestrial objects in orbit around WDs. 
In Sections 1.1, 1.2 we briefly discuss theoretical studies 
concerned with the likelihood of sub-stellar and planetary 
survival to stellar evolution. 



1.1 Sub-stellar companions to WDs 

Observationally, sub-stellar compani ons to WDs are fo und 
to be rare. Using the 2MASS survey, iFarihi et al l (|2005l ) es- 
timated that < 0.5% of WDs have L dwarf companions. 
More recently, excess near-infrared emission from WDs in 
the UKIDSS survey (Steele et al. in preparation) tentatively 
suggests the fraction of unresolved brown dwarf companions 
(including T dwarfs) may be slightly higher, between 1 — 2%. 
However, at the time of writing only three wide WD+BD 
systems have been spectro scopically confirmed, GD 16 5 
jBecklin fc Zuckermadfl98sl ) . PHL5038 dSteele et al.ll2009l ). 
and LSPM 1459 + 0857 AB (|Dav- Jones et all l20ld )~and 
two detached, non-eclipsing, short-period WD+BD sys- 
tems are currently known, WD0137 — 349 (|Maxted et alj 
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2004 iDobbie et al.ll2005l . Burleigh et 
al., in preparation, P « 9.9h). GD1400B and WD0137 - 
349B are the only two sub-stellar companions known to 
have survived the common envelope (CE) phase of stellar 
evolution, with WD0137 — 349B currently the lowest mass 
(~ 50M jup ) object known to have done so. 

Although infrared sky surveys such as UKIDSS, VISTA 
and WISE, and observatories such as Spitzer hope to reveal 
many more such binaries, they remain difficult to identify 
either as infra-red excesses or through radial velocity mea- 
surements. The detection of more close systems will allow us 



to place observational upper limits on the mass of sub-stellar 
companions that can survive CE evolution. Furthermore, 
examples of eclipsing WD+BD binaries will be important 
for e xploring the WD a nd substellar mass-radius relations 
(e.g. iParsons et afl|201Cl ). 

In addition the detection of a significant number of 
eclipsing WD+BD binary systems might help uncover the 
hypothesised population of 'old' cataclysmic variables (CVs) 
in which the current accretion rate is extremely low and 
the companion has been reduced to substellar mass (e.g . 
|PattersorJ[l99sl ;l Patterson et al ] |2005l ; lLittle"fair" et al 1 12003 ). 
While these systems elude direct detection as X-ray sources 
and remain difficult to identify in optical and infra-red sur- 
veys, it is possible to mea sure the mass and th e radius of 
the donor in eclipsing CVs. ILittlefair et all (|2006h confirmed 
the first such system through eclipse measurements, while 
ILittlefair et~ai1 |2007f ) showed that another eclipsing CV, 
SDSS J150722.30 + 523039.8, was formed directly from a 
detached WD/BD binary. Old CVs are important for shed- 
ding light on models of close binary evolution as well as for 
placing constraints on the period distribution of cataclysmic 
variabl es; in parti c ular, the period gap and t he period min- 
imum (|Kindll988l ; IParthasarathv et al.ll2007f l. 

1.2 Can planets survive stellar evolution? 

Every star less massive than 8 Mq (~ 97% of all stars in 
our galaxy) will end its life as a WD. Thus, it is natural 
to ask what will be the fate of known extra-solar planetary 
systems? This question also has particular interest for us, 



evolution is uncertain ( 
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cuss post-main sequence evolution of planetary systems 
and show that planetary survival is not beyond 
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locity observations of red giants indicate that planets or- 
bits beyond the radius of the star's envelope can sur- 
vive stellar evolut i on to that stage (see iFrink et al. 1 120021 ; 
lHatzes et al.ll2005l . ISato et al.ll2003l ). However, direct imag- 
ing searches at WDs have so far failed to detec t any plan- 
etary mass companio ns _(e.g. iHogan et all 120091 ). More re- 
cently, ISilvotti et al.l (|2007T l reported the detection of a ~ 
3 M jup planet orbiting an ex t reme horizontal branch star. 
Furthermore, iMullallv et all (|2008l ) found convincing evi- 
dence of a 2 M jup planet in a 4.5 year orbit around a pul- 
sating WD. The latter, if confirmed, will be the first planet 
detected in orbit around a WD, and will show that planets 
can indeed survive the death of their parent star. 

The existence of short-period planetary companions to 
WDs may seem less likely. Two scenarios may give rise to 
planets in short-period orbits around WDs: 
1) planets undergo CE evolution and survive their parent 
stars' evolution to a WD, or 2) their orbits are significantly 
changed by a process occurring at the end of the asymptotic 
gia nt branch (AGB) phase of stellar evolution. 

IVillaver fc Livid (|2007T ) investigated the fate of a planet 
engulfed by the envelope of an AGB star and suggested that 
planets in orbit within the reach of the AGB envelope will ei- 
ther totally evaporate or in rare cases, a more massive body 
may accrete mass and become a close companion to the star. 
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In this scenario only the massive companions (e.g. brown 
dwarfs like GD1400 and WD0137 - 349B) are likely to sur- 
vive the red giant and the AGB phases of stellar evolution. 
However, estimates of the minimum substellar mass neces- 
sary for survival are highly uncertain and depend on sev- 
eral f actors, for example, the efficiency of the envelope ejec- 
tion (|Villaver fc Liviq |2007| and references therein) . None 
the less, it is unlikely that terrestrial planets can surviv e 
engulfment and evaporation l|Wickramasinghe et al.ll201Cj ). 

Planets that escape engulfment by the red giant or 
asymptotic giant and that are sufficiently far from the stel- 
lar surface that they do not experience tidal drag, will have 
their orbital rad ii increased t o conserve angular momen tum 
(as described by |jeandll924r ).lD uncan fe Lissauerl (|l99St ) in- 
vestigated the stability of planetary systems during post 
main-sequence evolution, and found that for WD progen- 
itors experiencing substantial mass loss during the AGB 
phase, pla netary orbits become unst able on timescales of 
^10 year. iDebes fc Sigurdssonl (|2002f ) also studied the sta- 
bility of planetary systems and found that mass loss from 
the central star is sufficient to destabilise planetary sys- 
tems comprising two or more planets. For unstable systems 
in wh ich the orbits happen to cross, IDebes fc Sigurdssonl 
found that the most likely result was that one planet 
would be scattered into an inner orbit, while the other would 
either be boosted into a larger orbit, or ejected from the 
system altogether. This may result in WD systems which 
have settled into a configuration wherein planets are found 
at orbital radii which were originally occupied by the (now 
evaporated) inner planets before the RGB phase of stellar 
evolution. 

The above scenario provides a plausible explanation for 
the recent detection of silicate-rich dust discs around a grow- 
ing number of WDs at orbital radii up to ~ 1 Rq (e.g. 
Jurall2003l:|Reach et al.ll2005l : iFarihi et al]l20Qgl : IFarihi et all 
2009M Jural (|2003T ) suggests that the formation of dust discs 



around WDs is most probably due to the tidal disruption of 
an asteroid or larger body which has strayed too close to the 
parent star. Dynamical instabilities during the final stages 
of solar system evolution could have c aused the rocky body 
to m igrate inwards (as suggested by IDebes fc Sigurdssonl 
I2002T I. If the body wanders too close to the Roche radius 
of the WD it will be completely destro yed, produc ing a de- 
bris disc reminiscent of Saturn's rings l|jurall2003l). Recent 
studi es of the dust disc around the WD GD 362 ( Jura et all 
I2009T ) suggest that the more likely scenario which simultane- 
ously explains all of GD 362's distinctive properties is that 
we are witnessing the consequences of the tidal destruction 
of a single body that was as massive as Callisto or Mars. 

Consideration of dynamical interactions and orbital sta- 
bility indicates that while a terrestrial body may be per- 
turbed from some wide orbit into an eccentric orbit that 
takes it within the Roche radius of a WD and hence be 
disrupted, it is highly unlikely for such an object to be cap- 
tured into a stable, close orbit just beyond the Roche radius. 
None the less, one may speculate on the existence of 'shep- 
herd moons' accompanying the dust discs detected at WDs, 
similar to those at Saturn's rings. Alternatively, when close, 
double WDs are drawn together by gravitational radiation 
and merge, second generation terrestrial planets may form 
in remnant discs left by the tidal disruption of the lower 
mass degenerate jHanseiJliobl , iLivio et al.||2005l ). Indeed, 



IWickramasinghe et al.l |201Cf ) speculates that such an object 
may be closely orbiting the unusual magnetic WD GD 356. 
Thus, the existence of asteroids, moons and rocky planets 
in close orbits to WDs may not be entirely unreasonable. 

The detection of short-period substellar and planetary- 
mass companions to WDs, will open an exciting chapter in 
the study of extra-solar planet evolution, constraining the- 
oretical models of common envelope evolution and helping 
us to understand the ultimate fate of hot Jupiter systems as 
well as the fate of our own solar system in the post main- 
sequence phase. In this work we present the results of a 
study designed to investigate the detection limits for tran- 
siting sub— stellar and terrestrial companions in close orbits 
around WDsQ. In Section 2 we describe briefly the WASP 
project and WASP telescopes, which provide the WD light 
curves which form the basis of our transit search. In Sec- 
tion 3 we outline our Monte Carlo simulations, describing 
our detection method, as well as characterising the type of 
systems we might hope to detect. In Section 4 we discuss the 
results of our simulations and provide transit recovery rates 
for simulated light curves comprising random Gaussian noise 
(white noise) and correlated noise (red noise). In Section 5 
we present the results of a comprehensive transit search in a 
sample of 194 WD light curves found in the WASP archive. 
Finally, our conclusions are presented in Section 6. 



2 THE WASP PROJECT 

The WASP (Wide- Angle Search for Planets) project, oper- 
ates two robotic telescopes, one located amongst the Isaac 
Newton Group of telescopes (ING), in La Palma Spain, with 
a second instrument situated at the South African Astro- 
nomical Observatory (SAAO). Each instrument consists of 
eight f/1.8 Canon lenses each with an Andor CCD array of 
2048 2 13.5/j.m pixels, giving a field of view of 7.8 degrees 
square for each camera. The observation strategy is to cycli- 
cally raster the sky in a series of fields centred on the current 
local sidereal time and separated by 1 hour in right ascen- 
sion. Each observation lasts for about 1 minute (30 seconds 
exposure, plus slew and telescope settling time). This strat- 
egy yields well sampled light curves with a typical cadence 
of about 8 minutes per field. WASP provides good quality 
photometry with accuracy ^ 1% in the magnitude range 
V ~ 9-12. The WASP telesc opes and data analys is strate- 
gies are described in detail in lPollacco et al] (|2006r ). 



3 DETECTABILITY OF ECLIPSES AND 
TRANSITS OF WHITE DWARFS 

To assess the chances of detecting eclipses and transits of 
white dwarf host stars in ground-based wide-field surveys 
we performed an extensive set of Monte Carlo simulations. 
The approach we adopted was to create realistic synthetic 
light curves containing eclipse and transit signatures of the 
expected depth and duration for a range of companion size 
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of discovering asteroids and moons in much wider orbits around 
WDs in Kepler data. 
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and orbital period, then to attempt to detect these sig- 
natures using a standard transit-detection algorithm (box- 
least-squares, BLS). By noting the rate at which the BLS 
search recovered the transit at the correct period (or an 
integer multiple or fraction) we were able to estimate the 
feasibility of detecting such systems in an automated man- 
ner. 



3.1 Characteristics of the transit signal 

The probability (ptr) that a low mass star, brown dwarf or 
planet in a circular orbit will transit or eclipse its host star 
is given by 
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Assuming an orbital inclination i — 90° the depth (St, 
and duration (D tI ) of such a transit are given by 
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For our simulations we have chosen the parameters of 
the host star to represent a typical 1- Gyr-old carbon-core 
WD of mass Af» = 0.6 M Q and radius R, = 0.013 R . We 
explored the detectability of planetary transits across the 
two-dimensional parameter space defined by the orbital pe- 
riod and the planet radius. We considered orbital periods in 
the range P ~ 2 hours to 15 days (equivalent to orbital dis- 
tances of between a ~ 0.003 and 0.1 AU). The lower-limit to 
the orbital period was chosen to yield an orbital separation 
close to the Roche radius of the WD, the upper-limit by a 
requirement that we have a reasonable chance of detecting 
five or more transits in a typical 150 day observing season 
of a WASP survey field. 

Figure [1] shows the probability that a given system will 
transit, and the depth and duration of such transits across 
this parameter space. It is evident from this diagram that the 
signatures of transits of WDs by typical planet-sized bodies 
will be rather different than those seen for typical transit- 
ing hot Jupiters. In particular the transit duration is much 
shorter for WDs than for normal stars (from ~ 1-30 min for 
companions with sizes ranging from Moon-size to Jupiter- 
size, compared to 2-3 hours for a typical hot Jupiter), and 
the transit depths are much larger (from around 3% for 
a Moon-sized to 100% for any companion larger than the 
Earth, compared to ~ 1% for a hot Jupiter). 



3.2 Generation of synthetic light curves 

The synthetic light curves were generated using the time 
sampling of a typical WASP survey field, light curves were 
generated with statistical signal-to-noise ratios representa- 
tive of three magnitude ranges (V ~ 12, 13 and 15) span- 
ning the range of brightness of WDs in the WASP survey. 
The corresponding photometric accuracy of WASP over this 
range is ~1% to 10%. Because WASP data show resid- 
ual covariant-noise structure due to instrumental system- 
atics we have tested the transit recovery rate in the case 
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Figure 1. Contours of constant transit probability (top), depth 
(centre) and duration (bottom) in the parameter space defined by 
orbital period and planetary radius. The transit probability and 
depth are expressed in percentage values. The transit duration is 
expressed in minutes. 



of both uncorrelated "white" noise and correlated "red" 
noise. In the white-noise case we injected transit-like sig- 
natures into otherwise non- variable light curves, adding a 
Gaussian-distributed noise component of standard devia- 
tion a. We chose a to be representative of the mean pho- 
tometric error on the points obtained from a real WASP 
light curve for an object of our chosen magnitude. The tem- 
plate WASP light curves therefore defined the time-sampling 
and the average signal-to-noise, but the photometry was 
otherwise entirely synthetic. In the red noise case we in- 
jected fake transits into a set of unmodified WASP light 
curves obtained from a densely sampled field observed dur- 
ing the 2004 season. Data from the 2004 season have been 
detrended and t horoughly searched for transi t-like events 
as described in ICollier Cameron et al.l (|2006T l. Moreover, 
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Figure 2. Two example synthetic light curves. Top, an eclipsing 
brown dwarf in an orbit with a period of 2hr, around a WD. 
Bottom, a 1.2 companion to a WD in 5hr orbit. 



we cross-correlated stars in the WASP field with the pub- 
licly available Gen eral Catalogue of Variable Stars (GCVS; 
ISamus et al.ll20o3 ). and we evaluated the RMS of each light 
curve which we used to identify and remove variable objects 
after individual eye-balling. Finally, each light curve con- 
tained around 4240 data points, acquired over 116 nights, 
and spread across a baseline of 128 nights, with photomet- 
ric accuracy raging from ~ 1% to ~ 10% for stars in the 
magnitude range 12 < V < 15 . For each light curve we 
used the WASP pipeline fluxes and errors derived after de- 
trend ing by the SysRem algorithm (jCollier Cameron et al.l 
l200rj ). 

Planet transit light curves of main-sequence stars show 
a characteristic shape, with an ingress lasting several tens of 
minutes, a flat bottom of 2-3 hours and an egress again last- 
ing tens of minutes. For the case of a WD host star consid- 
ered here, the ingress and egress duration is typically short 
compared to cadence of the WASP survey (8-10 minutes). 
We therefore ignore the detailed shape of the ingress and 
egress phases and modelled the transit signatures as simple 
box-like profiles. 

To cover the orbital period-planet radius parameter 
space we selected seven trial periods spaced approximately 
logarithmically (P = 0.08, 0.22, 0.87, 1.56, 3.57, 8.30 and 
14.72 days), and five planet radii R p = 10.0, 1.0, 0.6, 
0.34 and 0.27 R®. We modelled the set of synthetic light 
curves by injecting fake transit signals into phase-folded 
light curves at the trial period with a random transit epoch 
to in the range < to < P. We computed the transit du- 
ration according to Equation O and hence the width of the 
transit in orbital phase cV = Dt r /P. For all data points 
falling in the phase range ^ (j>i ^ cj>t T we then reduced 
the observed flux by a factor <5 tr . For each combination of 
orbital period and planet radius we generated 100 synthetic 
light curves. 

Figure [2] shows two examples of our simulated transit 
light curves. The top panel shows the synthetic light curve 
of an hypothetical eclipsing WD+BD binary system with 
an orbital period of P = 1 16 mins, similar to W D0137 — 349 
(a non-eclipsing system, iMaxted et al.l 120061 ) . The lower 



panel shows the simulated transit light curve for a rocky 
body of radius 1.2 R® in a 5hr orbit. 



3.3 Detection algorithm 

To recover the transit signals from the synthetic light curves 
we used an implementation of the box-least-squares (BLS) 
algorithm (|Kovacs et al.l |2002| ) commonly used to detect 
transits of main sequence stars. The BLS algorithm is most 
sensitive when the modelled box-width closely matches the 
duration of the true transit signal. Thus, to ensure that 
the BLS search was sensitive across the expected range of 
transit durations, we chose to search a grid of box widths 
Wb = {1, 2, 4, 8, 16, 32} minutes, covering the range in tran- 
sit durations over most of our parameter space (Figure [T|. 
We defined the grid of trial periods sampled by BLS as fol- 
lows: 



iv 



(4) 



where P m i n = 2hr and P max = 15 day. The frequency inter- 
val was chosen such that the accumulated phase difference 
between successive trial frequencies over the duration of the 
light curve corresponds to the width of the shortest trial 
box duration at the longest period searched. At each trial 
frequency we defined a set of trial transit epochs at an inter- 
val chosen such that W{, — Wb for the shortest trial box 
duration, adjusted to meet a constraint that the number of 

epochs N e = P/W^ be an int eger. 

Adopting the notation of ICollier Cameron et all ((2006) 
we denote the set of observations in the light curve £i with 
formal variances erf and additional variances crf^ computed 
by SysRem to account for transient systematic variations 
due to patchy atmospheric extinction, for example. For each 
data point we compute a weight 



(i) 

then subtract the weighted mean of the observations 
J2i £iWi 



t = 



to obtain Xi — x~i 



x. We then define 



summing across the whole dataset. At each trial period 
we fold the light curve and accumulate into a set of bins j 
of width W^ the following quantities 



Sj — ^ ^ XilOif 

The fitted transit depth in the bin, and its associated 
variance are then 

Var(<V 



and the signal-to-noise ratio of a putative transit in the 
bin is: 



<S, = 



V /Var^) 
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Figure 3. BLS power spectrum for the transit signal of an Earth- 
sized body in an orbit with a 5 hour period. While the correct 
(inserted) period is recovered (as indicated by the dashed line) by 
the standard implementation of the BLS algorithm, strong alias- 
ing and noise structure are present, particularly at long periods. 



For each epoch bin we compute the signal-to-noise Sjfi 
for each of our trial box widths by co-adding values of Sj 
and Tj in adjacent bins and find the maximum Sj, ma x over 
all of the box widths. The maximum value of Sj, ma x across 
all of the epoch bins then represents the significance of the 
detection of a transit at the given trial period. When com- 
puting this latter maximum we ignored trial boxes which 
contained less than five data points, or did not contain data 
points from at least five distinct orbits. 

Figure [3] shows the periodogram computed in this fash- 
ion for the sample light curve shown in the lower panel of 
Figure [2] a 1.2 R® radius body in a 5 hour orbit. Although 
the correct simulated period shows up as the highest peak 
in the periodogram, there is significant non-random struc- 
ture in the noise continuum, particularly for trial periods 
longer than 1 day. Our interpretation of this phenomenon is 
that it is a by-product of a few unique features of the tran- 
sit signals we are dealing with. Firstly the duration of the 
transits, and hence the width of the boxes fitted by the BLS 
algorithm, is much shorter with respect to both the orbital 
period and the WASP survey cadence than for transits of 
main sequence stars. As a consequence the trial bins will 
contain many fewer data points than for a main sequence 
transit search, particularly at longer trial periods. Secondly 
as the transit signals are so deep compared to the main se- 
quence case, they are less prone to being "washed out" when 
the light curve is folded on an incorrect trial period and the 
in-transit points spread across all orbital phases. In the case 
illustrated the Earth-sized body will cause transits with a 
depth of ~ 70%. Even in a trial bin with, for example, 1 
in-transit and 10 out-of-transit data points the presence of 
the single in-transit point would drag down the mean light 
level in the bin by more than 6%, which could be sufficient 
to be regarded as a significant detection. 

To address this issue we modified equation [5] as follows 

S'i = . 5j (6) 
JVar(5 i ) + ^ 



Figure 4. BLS power spectrum as in Figure l3l obtained with the 
improved BLS routine. For the same transit signal we achieve 
higher S/N values and higher statistical significance for the de- 
tection. The strong aliasing and the noise structure seen in Figure 
[3] is much reduced. 

where dj is the mean-square deviation of the data points 
within the bin about their mean. This modification will have 
the effect of strongly reducing the computed signal-to-noise 
of bins which contain a mix of in-transit and out-of-transit 
data points. The magnitude of this down-weighting will tend 
to increase as the depth of the transit signal increases, but 
only in the cases in which the trial period does not match 
the true period (or an integer multiple or fraction). Where 
a bin contains just in-transit points the down-weighting will 
generally be small for all transit depths. Figure 0] shows the 
BLS periodogram computed using this modified prescription 
for the signal-to-noise, for the same transit as in Figure [3] 

3.4 False detection rate 

Automated searches for weak signals in noisy data are inher- 
ently susceptible to "false alarms" , whereby a chance align- 
ment of noise fluctuations in the data are misinterpreted by 
the search algorithm as a evidence likely detection of the sig- 
nal being hunted. It is useful therefore to be able to define a 
filter which can be applied automatically to weed out these 
false detections. We achieved this by constructing synthetic 
light curves based on the time-sampling of sample WASP 
light curves which contain pure white-noise, but no simu- 
lated transit signal. We then computed BLS periodograms 
for these light curves in the same manner as for those con- 
tain ing simulated transi ts. 

iKovacs et al.l l| 2002l) define a useful metric for assessing 
the likely significance of a peak in a BLS periodogram, which 
they refer to as the Signal Detection Efficiency (SDE): 

SDE = 5pcak ~ - 

where <S pca k is the height of the peak, and S and 175 
are measures of the mean level and scatter in the noise 
continuum of the periodogram. We computed the SDE of 
the highest peak in each of the synthetic, transitless light 
curves. The cumulative distribution function of these SDE 
measures over the whole sample is plotted in Figure [5] This 
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X 

Figure 5. Probability Distribution Function (PDF) for the 
modificd-optimiscd BLS routine. The dashed-line shows the de- 
tection threshold 6.3SDE for a 10% noise contribution. 



distribution function allowed us to define a threshold SDE 
below which we could automatically discount a detection 
as likely to be a false alarm. As the size of the sample of 
WDs observed by WASP is fairly small, just a few hundred, 
we chose a relatively generous threshold (SDE t hrcsh = 6.3) 
which would allow through around 10% of false detections. 
For larger surveys a more strict threshold might be necessary 
to avoid being swamped by false detections. 



3.5 Recovery rates of synthetic transits 

Figure [6] and Tables [1] [2] and [3] summarise our recovery 
rate for simulated transit signals injected into synthetic light 
curves of WDs of magnitudes V ~ 12, V ~ 13 and V ~ 15 
respectively. We regard as a match any trial in which the 
most significant detected period is within 1% of being an 
integer fraction or multiple from 1/5 x to 5x the injected 
transit signal. 

We have attempted to separate out the various factors 
which can affect the efficiency of detection of these transit 
signals. When generating each synthetic light curve we can 
readily assess a priori whether it will fail the tests requiring 
a minimum number of individual transits and in-transit data 
points. We list in Tables \T\ [2] and the fraction /flit which 
pass these two tests. It is evident from these tables that these 
requirements alone render transiting companions essentially 
undetectable at our longest trial periods (8.30 and 14.72 
days) in a WASP-like survey; the transits are too short in 
duration and too infrequent to be adequately sampled. For 
companions around 1 R® and larger however there is a good 
chance of detection out to periods of around 4 days, at least 
in principle. 

The table also shows the impact of adding representa- 
tive photometric noise on the detection rates (/det)- For the 
idealised photon-noise-limited case objects as small as Mer- 
cury could be detected to periods of around 1.5 d, and the 
Moon for periods less than 1 d. Once the impact correlated 
instrumental noise (red noise) is added, Moon-sized com- 
panions become almost undetectable, though the recovery 
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Figure 6. Recovery rate of simulated transits of a WD of mag- 
nitude V ~ 12 (top), V ~ 13 (middle) and V ~ 15 (bottom). The 
frequency contours are expressed in percentage values. 



rates for larger bodies, particularly in short-period orbits, 
remains encouraging. 

Our key conclusion from these simulations is that for the 
case of transits of WDs, the degree of photometric precision 
delivered by a survey is of somewhat secondary importance 
compared to a high cadence and continuous coverage. For 
planet-sized bodies individual transits will be quite deep and 
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Table 1. Recovery rate of simulated transits of a bright WD (V ~ 12). Results are shown for synthetic light curves containing white 
and red noise. In both cases f^ et is the fraction of cases in which the highest peak in the pcriodogram satisfies our period matching 
criteria and has an SDE > 6.3, and /t> t is the fraction which match the period criteria but have SDE < 6.3. /fli t is the fraction of cases 
in which the synthetiscd light curves pass the requirements for a minimum number of transits (> 5) and data points (> 5) in transit. 
Dashes indicate cases in which the companion would be tidally disrupted within the WD's Roche radius. 
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readily detectable in data of moderate photometric quality, 
however it is the short duration of the transits that is the 
main factor limiting the transit detection rate in surveys 
optimised for main sequence stars. 



4 SEARCHING FOR TRANSIT SIGNALS IN 
WASP SURVEY DATA 

Encouraged by the results of our simulations we selected a 
sample of WDs, which have been routinely monitored by 



WASP through the 2004-2008 observing seasons, and per- 
formed a systematic search for eclipsing and transiting sub- 
stellar and planetary companions. We selected the sample by 
cross-correlating the catalogue of WASP objects for which 
more than 600 data points are availab le with the McCook & 
a ion catalogue jMcCook fc Sionl2003l ). The resulting sample 
of 194 WDs with magnitude V < 15 is presented in Table 4. 

We searched the sample for transits and eclipses using 
our implementation of the BLS algorithm, searching peri- 
ods ranging from 2 hours to 15 days. In addition we have 
also inspected each of the individual light curves by eye. In 
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both searches we found no evidence for any transiting and 
eclipsing companions within the period range searched in 
this study. We have used this null result together with the 
results of our simulations to estimate an upper-limit to the 
frequency of such close companions for the sample of WDs 
considered in this study. 



4.1 



Limits on frequency of companions to WDs 



In order to estimate an upper limit to the frequency of close 
substellar and planetary companions to WDs, we used the 
detection limits derived from our simulations and the re- 
sults obtained from the analysis of the sample of 194 WDs. 
We first used a binomial distribution to describe the prob- 
ability V(n\ N, /) of finding n transiting companions for a 
given sam ple of TV stars, with a true com panion frequency 
f (e.g.. see McCarthy fc Zuckermanll2004l ; and Appendix of 
iBurgasser et alj 20031 ') as follows: 

N\ 



V(n;N,f) = 



nl(N-n)l 



f (1 ^ fY 



(7) 



When the two quantities TV and n are known equation [7J 
can be used to derive the distribution (Pi) describing the 
probability of /, where / is the frequency of transiting 
companions. The probability Vi(f;n,N) is proportional to 
V(n; N, /) for / in the interval [0, 1]. We obtain V\ by nor- 
malising : 



/ 



Pi(f;n,N) df 



(8) 



which yields Vi = (N + l)V. 

Although our complete sample numbers N = 194 stars, 
we have already established that even if all of these have 
companions only a fraction p t r(R P ,P) will exhibit a tran- 
sit, and of those which do exhibit a transit only a frac- 
tion pdct(R P , P) would be detectable in a WASP-like survey. 
Both of these factors will act to reduce the total number of 
transiting companions detected in the survey, or in the case 
of a null result will tend to weaken the constraints that can 
be placed on true companion frequency by such a survey. 
To incorporate these factors we modified our effective sam- 
ple size as: 

N' = N x P tr(R P , P) x Pdct (R P , P) 

and used this in Eq. [51 which we integrated to find 
the limiting companion frequency /ii m that encloses 95% 
of the probability distribution. Figure [JJ (top panel) shows 
the upper-limit on companion frequency for a null detection 
in a "perfect" survey in which pact = 1 and a sample size 
N = 194. In such a survey the detectability of companions is 
limited solely by the intrinsic probability of them transiting 
their host. 

To factor in the efficiency of detection of transits in 
a WASP-like survey we need to determine a representative 
Pdet{Rp, P)- Our simulations were performed at only three 
specific host-star magnitudes, whereas the distribution of 
the magnitudes of the stars in our sample is of course a 
continuum (covering the range V ~ 9 — 15). We therefore 
combine the three magnitude-specific pdct maps (Figure [6]) 
into a single map by interpolating/extrapolating according 
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Figure 7. Top panel: upper-limit on companion frequency in- 
ferred from a null detection in a survey sample of size N = 194 
assuming perfect detectability of transits across the parameter 
space. Lower panel: limits in the same sample folding in the de- 
tectability of transiting systems in a WASP-like survey. In both 
panels the frequency contours are expressed in percentage values. 



to the magnitude of each object in our sample and com- 
bining these to form an averaged map which can be folded 
in to our calculation of the upper-limits. The resulting lim- 
its corresponding to the 95% of the integrated probability, 
are shown in in the lower panel of Figure [Jj Our results 
show that for rocky bodies smaller than the size of Mer- 
cury no useful upper limits to the frequency of companions 
to WDs can be found, and that for Earth-sized companions 
only weak constraints can be imposed. However, it does sug- 
gest that objects the size of brown dwarfs or gas giants with 
orbital periods P< 0.1 — 0.2 days must be relatively rare 
(upper limit of ~10%). 



5 CONCLUSION 

We have used a modified version of the Box Least Squares 
algorithm to investigate the detection limits for substellar 
and planetary companions to WDs achievable using data al- 
ready available in the WASP photometric survey. Our simu- 
lations proved extremely encouraging, suggesting that plan- 
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etary bodies as small as Mercury at small orbital radii can 
be detected with good photometric data even in the pres- 
ence of red noise. For smaller bodies red noise in the light 
curves becomes increasingly problematic, while for bodies 
with larger orbital periods, the absence of significant num- 
bers of in-transit points, significantly decreases the detection 
sensitivity. 

Application of our modified BLS algorithm to search for 
companions to WDs in a sample of 194 stars in the magni- 
tude range V ~ 9-15, available in the WASP archive, did 
not reveal any eclipsing or transiting substellar or planetary 
companions. Visual inspection of individual light curves for 
the WDs in our sample confirmed the absence of significant 
periodic dropouts in the WASP data. We have used the non- 
detection of planetary companions to the WDs in our sam- 
ple together with the estimated detection sensitivities de- 
termined from our simulations, to place upper limits to the 
frequency of substellar and planetary companions to white 
dwarfs. While no useful limits can be placed on the likely 
frequency of Mercury-sized or smaller companions, and only 
weak constraints on the frequency of Earth-sized objects in 
the closest orbits, slightly stronger constraints can be placed 
on the frequency of larger bodies in very short-period orbits. 
For example, brown dwarfs and gas giants radius ~ R jup 
with periods < 0.1 — 0.2 days, similar to the known WD+BD 
binary WD0137 - 349 (jMaxted et al.ll200rj ). must certainly 
be relatively rare (< 10%). Of course, this limit needs to be 
compared with those derived from other s o urces , e.g. infra- 
red sky surveys. For example, iFarihi et alj (|2005h estimated 
that < 0.5% of WDs have L dwarf companions, while Steele 
et al. (in prep.) tentatively suggests the fraction of unre- 
solved brown dwarf companions (including T dwarfs) may 
be slightly higher, bet ween 1 — 2%. From Spitzer photome- 
try IFarihi et ail (120081 ) suggests that < 4% of white dwarfs 
have unresolved substellar companions > 10M jup , although 
the limits at lower masses (e.g. < 6M jup ) are considerably 
weak. 

To place more stringent constraints on close substel- 
lar and gas giant companions to WDs, and similarly strin- 
gent constraints on Earth size bodies in close orbital sep- 
arations likely requires significantly larger WD samples. In 
addition, our simulations and analysis of WD light curves in 
the WASP archive suggests the degree of photometric pre- 
cision is of somewhat secondary importance compared to a 
high cadence and continuous coverage. The short duration 
of eclipses and transits of WDs (~ 5 — 20 mins for com- 
panions radius ~ Rj up ; ~ 1 — 5 mins for terrestrial bodies) 
compared to the ~ 8 min cadence of WASP observations, 
appears to be the main factor limiting the transit detection 
rate in a survey optimised for planetary transits of main 
sequence stars. 

Future surveys such as Pan-STARRS and LSST will 
be capable of detecting tens of thousands of WDs. How- 
ever, we emphasise that observations of high cadence and 
long baseline are of greatest benefit when attempting to de- 
tect the signature of close, eclipsing and transiting substellar 
and planetary comp anions to WDs. Space missions such as 
COROT, Kevler (see bi Stefano et aljboid ) and, especially, 
PLATO may therefore be better suited to a survey of WDs 
as they deliver uninterrupted coverage at high cadence and 
exquisite photometric precision (~ 10 -4 — 10 -5 ) and could 



at least in principle detect the transits of asteroid-sized bod- 
ies across a WD. 
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Table 2. Recovery rate of simulated transits of a V ~ 13 WD 



White noise red noise 



Size Rpi 5ti P Dtr / fllt / dct f ht / dct f ht 

R (%) (days) (min) (%) (%) (%) (%) (%) 
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Earth 1.0 49 



Mercury 0.45 10 
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Table 3. Recovery rate of simulated transits of a V ~ 15 WD 



White noise red noise 



Size R p i 5ti P Dtr / fllt / dct f ht / dct f ht 

R (%) (days) (min) (%) (%) (%) (%) (%) 



BD-Gas Giant 10.0 100 



Earth 1.0 49 
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Moon 0.27 3.6 
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WDs observed by WASP, including the WASP identity, corresponding identity in the McCook & Sion catalogue, WASP 
magnitude, and the number of individual data points contributing to the light curve in the WASP archive. The WASP 
magnitude is defined as-2.51ogio(-F/10 6 ), where F is the mean WASP flux in //Vega; it is a pseudo-V magnitude comparable 

to the Tycho-V magnitude. 
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